Modified-Puls  Routing 
in 

Chuquatonchee  Creek 


DTIC 

ELECTE 
MAY  3  9  1985 


Research  Document  No.  23 


September  1980 


This  document  has  been  approved 
fcr  public  release  and  rale;  its 
d'  t  ib-j’ien  is  unlimited. 


F 


85  5  23  12i 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dele  Entered) 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 

Research  Document  No.  23 


4.  TITLE  (end  Subtitle) 

Modified-Puls  Routing  in  Chuquatonchee  Creek 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 

s' 


5.  TYPE  OF  REPORT  &  PERIOD  COVERED 


S-  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORS 


8.  CONTRACT  OR  GRANT  NUMBERf.J 


Theodor  Strelkoff 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


4.  MONITORING  AGENCY  NAME  a  ADDRESS (tt  dllterent  from  Controlling  Of/Ice) 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  a  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 


April  1985 


13.  NUMBER  OF  PAGES 


ISO.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thle  Report) 


Distribution  of  this  document  is  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abatract  entered  In  Block  20,  It  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 

This  work  was  performed  under  contract  (DACW05-80-P-0324)  to  the 
Hydrologic  Engineering  Center. 


19.  KEY  WORDS  (Continue  on  rarer  aa  aide  If  necemeary  and  Identity  by  block  number) 

Flood  Routing,  Modified-Puls,  Storage  Routing,  Hydrologic  Routing 


20.  ABSTRACT  (Continue  an  rarer* 


r  ft  nejcaaaeey  and.  Identity  by  block  number) 


The  Modified-Puls  method  has  at  its  core  the  postulate  that  storage 
depends  only  on  outflow  rate.  The  relationship  S(0)  is  preseined  known  for  the 
reach  of  interest  and  for  the  range  of  outflow  to  be  encountered.  The  hope 
that  the  method  will  yield  at  least  approximately  correct  outflow  hydrographs 
rests  on  the  assumption  that  the  storage  depends  primarily,  if  not  only,  on 
outflow  rate.  If  only  certain  features  of  the  outflow  hydrograph,  such  as 
peak  and  time  of  peak,  are  desired  with  reasonable  precision,  then  the 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfHTiwi  Palm  Enfrfd) _ 

dependence  of  S  on  0  perhaps  could  be  relaxed  further,  with  dependence 
existing  only  for  some  particularly  important  range  of  0. 

>- 

The  given  reach  of  the  Chuquatonchee  was  schematicized  by  a  simple 
prismatic  geometry,  with  a  rectangular  main  channel  for  in-bank  flows,  to 
receive  the  large  flood  flows.  The  unsteady  flows  generated  by  the  test  floods 
in  the  test  channel  were  modeled  by  solution  of  the  complete  one-dimensional 
Saint-Venant  equations.  Water-surface  profiles  (and,  hence,  stored  volumes), 
as  well  as  outflow,  were  monitored  during  the  course  of  the  calculation.  Time 
and  distance  steps  in  the  numerical  procedures  were  taken  small  enough  that 
their  size  played  no  significant  role  in  the  resulting  solutions. 

The  modified-puls  method  with  a  storage-outflow  based  on  steady-flow 
profiles  is  not  suitable  for  this  reach,  primarily  because  of  the  great  size 
of  the  depth  gradient  relative  to  other  forces  driving  or  retarding  the  flow, 
and  because  the  depth  gradient  is  strongly  a  function  of  the  unsteadiness  in  a 
channel  of  small  slope.  Further,  the  effect  of  depth  on  storage  is  highly 
exaggerated  by  the  broad  flat  floodplains  of  this  reach.  By  the  same  token, 
S(O)  is  highly  event  dependent,  so  that  no  single  function,  even  approximately 
valid  for  a  wide  range  of  peaks,  exists. 

While  doubt  has  been  casL  on  the  accuracy  of  the  topographical  data 
gathered  for  Chuquatonchee  Creek,  it  was  found  that  the  errors  stemming  from 
application  of  Modified-Puls  to  events  in  reaches  for  which  it  is  theoretically 
unsuitable,  were  of  sufficient  magnitude  that  they  alone  could  explain  the 
noted  discrepancy. 
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PREFACE 


The  work  reported  herein  is  a  product  of  the  Hydrologic  Engineering  Center's 
continuing  effort  to  improve  both  the  accuracy  and  efficiency  of  the  analytic 
techniques  used  by  the  Corps  of  Engineers  to  route  floods  through  natural  and 
modified  river  channels.  This  report  describes  the  application  of  a  technique 
for  evaluating  the  relative  accuracy  of  flood  routing  methods  to  a  natural 
river.  The  theoretical  development  of  that  technique  is  presented  in  the 
companion  report  "Comparative  Analysis  of  Flood  Routing  Methods." 
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1.  INTRODUCTION  AND  SUMMARY 


The  principle  of  mass  conservation  can  be  expressed  for  a  reach  of 
river  by  the  equation 

1  +  T  -  0  =  (1) 

in  which  the  terms  on  the  left  represent  volumetric  rates  of  flow,  I  the 
inflow  into  the  reach  at  its  upper  end,  T  the  instantaneous  sum  of  all 
tributary  inflow  over  the  length  of  the  reach,  *  and  0  the  outflow  at  the 
lower  end.  The  term  S  represents  the  volume  of  water  in  the  reach,  and  t  is 
time.  Equation  1  is  exact,  and  applies  at  every  instant  during  the  total 
time  of  interest.  With  the  inflow  hydrograph  I(t)  known,  if  also  T(t)  and 
S(t)  were  known,  determination  of  0(t) ,  the  outflow  hydrograph,  would  be 
trivial.  Of  course,  S(t)  for  any  given  inflow  hydrograph  is  not  known.  So, 
the  factors  influencing  S  and  a  mathematical  expression  of  the  dependency 
are  sought  to  allow,  ultimately,  solution  of  Eq.  1. 

The  Modif ied-Puls  method  has  at  its  core  the  postulate  that  storage 
depends  only  on  outflow  rate.  The  relationship  S(0)  is  presumed  known  for 
the  reach  of  interest  and  for  the  range  of  outflow  to  be  encountered.  Then, 
Eq.  1  takes  the  form 


if +  0-I  +  T 


(2) 


an  ordinary  differential  equation  for  the  outflow  hydrograph.  With  I(t) 

J  c 

and  -tI-(O)  known,  and  the  tributary  hydrograph  known,  assumed,  or  neglected, 
dO 

the  outflow  hydrograph  0(t)  can  be  computed.  The  Modif ied-Puls  method  pro¬ 
vides  a  numerical  technique  for  solving  this  differential  equation. 


^Understood  as  net  inflow:  tributary  inflow  +  overland  inflow  +  precipitation 
-  evaporation  -  seepage  -  lateral  outflow  through  side  weirs,  pump  intakes, 
diversions,  etc. 


The  hope  that  the  method  will  yield  at  least  approximately  correct 
outflow  hydrographs  rests  on  the  assumption  that  the  storage  depends  primarily 
if  not  only,  on  outflow  rate.  If  only  certain  features  of  the  outflow  hydro¬ 
graph,  such  as  peak  and  time  of  peak,  are  desired  with  reasonable  precision, 
then  the  dependence  of  S  on  0  perhaps  could  be  relaxed  further,  with  depen¬ 
dence  existing  only  for  some  particularly  important  range  of  0. 

A  key  factor  in  the  application  of  the  Modif ied-Puls  method  to  any 
given  reach  is  the  determination  of  this  relationship  S(0)  for  the  reach. 

Two  techniques  of  arriving  at  this  function  are  employed  extensively: 

(a)  computing  steady-flow  profiles  in  the  given  reach  for  a  succession  of 
discharges  0,  from  given  topographical  and  roughness  information,  then 
computing  water  volume  under  the  profile  for  each  value  of  discharge,  and 

(b)  deriving  S(0)  for  historical  events  from  given  inflow  and  outflow  hydro¬ 
graphs.  The  former  procedure  is  complicated  by  the  difficulty  of  obtaining 
accurate  topographical  and  roughness  data,  and  by  the  fact  that  these  can 
vary  from  event  to  event.  The  latter  is  complicated  by  the  existence  of 
ungaged  tributary  inflow,  with  only  /T  dt  derivable  from  the  inflow  and 
outflow  hydrographs. 


/ 

t 


T  dt  =  /  0  dt  -  /  I  dt 


(3) 


with  tj  and  t ^  taken,  at  the  beginning  and  end  of  the  event,  such  that  S(tj) 

=  S(t2)-  With  T(t)  and  S(0)  independent  functions,  Eq.  2,  alone,  is  insuf¬ 
ficient  to  solve  for  either,  if  only  I(t)  and  0(t)  are  known.  In  one  variant 
(reference  1),  an  "optimum"  S(0),  common  to  several  (three)  calibration 
events,  was  found  by  minimizing  the  negative  portions  of  the  resultant 
computed  T(t). 

Whichever  procedure  is  used,  it  relies  on  the  validity  of  the  assumption 


that  storage  is  truly  a  function  of  outflow  alone,  at  least  approximately. 


-3- 


In  the  given  reach  of  Chuquatonchee  Creek,  the  two  procedures  led  to 
very  different  S(0)  functions.  It  was  the  purpose  of  the  present  investi¬ 
gation  to  determine  the  most  likely  causes  of  the  disagreement. 

The  given  reach  of  the  Chuquatonchee  was  schematicized  by  a  simple 
prismatic  geometry,  with  a  rectangular  main  channel  for  in-bank  flows, 
and  a  wide,  plane,  slightly  sloping  floodplain  bordered  by  vertical  bluffs 
to  receive  the  large  flood  flows.  Specific  flood  events  of  record  were 
also  schematicized,  by  simple,  single-peak  inflow  hydrographs,  with  con¬ 
tinuously  variable  time-rate  of  change.  Test-channel  geometry  was  selected 
consistent  with  the  topographical  features  of  the  Chuquatonchee  valley. 

Peak  and  rise  time  of  test  inflow  hydrographs  were  consistent  with  peaks 
and  rise  times  of  selected  historical  events. 

The  unsteady  flows  generated  by  the  test  floods  in  the  test  channel 
were  modeled  by  solution  of  the  complete  one -dimensional  Saint-Venant 
equations.  Water-surface  profiles  (and,  hence,  stored  volumes),  as  well 
as  outflow,  were  monitored  during  the  course  of  the  calculation.  Time  and 
distance  steps  in  the  numerical  procedures  were  taken  small  enough  that 
their  size  played  no  significant  role  in  the  resulting  solutions.  Because 
test-channel  geometry  and  roughness  are  known  precisely,  errors  in  input 
data  are  not  of  concern,  and  solutions  of  the  Saint-Venant  equations  can 
be  considered,  except  for  one  reservation,  very  close  to  what  would  be 
observed  in  real  water  flows.  To  be  kept  in  mind  is  that  the  model  assumes 
lateral  flooding  of  the  floodplain  to  take  place  instantaneously,  in  keeping 
with  the  one-dimensional  character  of  the  model,  which  assumes,  at  all 
times,  a  horizontal  water  surface  transverse  to  the  main  channel.  In  actual 
fact,  water  flows  up  onto  the  floodplain  in  consequence  of  a  transverse 


sloping  water  surface,  and  a  finite  length  of  time  is  required  to  inundate 
the  plain  to  any  given  elevation. 
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In  an  opposing  set  of  calculations,  stored  volumes  were  determined  as 
a  function  of  outflow  for  a  succession  of  steady  states  at  various  discharges. 
The  resultant  S(0)  function  was  placed  into  a  numerical  solution  procedure 
for  Eq.  2,  and  Modified-Puls-routed  hydrographs  determined.  As  a  variation 
of  this  procedure,  routing  could  be  performed  over  a  series  of  subreaches, 
the  outflow  hydrograph  from  one  subreach  serving  as  the  inflow  hydrograph 
for  the  next  subreach  downstream.  The  S(0)  relation  for  each  subreach  was 
determined  by  apportioning  the  storage  determined  for  the  entire  reach 
amongst  the  subreaches  in  proportion  to  their  length.  For  the  uniform  flows 
developed  for  each  steady-state  discharge  (in  consequence  of  the  prismatic 
channel  and  normal-depth  rating  curve  postulated  at  the  downstream  end) , 
this  is  a  precisely  correct  distribution  of  storage. 

Tests  were  performed  without  tributary  inflow,  with  tributary  inflow 
distributed  uniformly  and  nonuniformly  over  the  length  of  the  reach,  and 
also  lumped  at  the  upstream  and  downstream  ends. 

The  most  significant  finding  was  that,  in  this  reach,  storage  is  only 
very  loosely  tied  to  outflow.  Mathematically  S  ^  S(0),  but  instead,  S  = 
S(0,....),  and  that  0  is  not  even  a  very  important  argument  during  much 
of  the  event.  Physically,  the  shape  of  the  flood  wave,  and  hence  the  storage 
therein,  is  very  different  from  the  shape  of  the  profile  in  steady  flow. 

The  key  feature  of  the  Chuquatonchee  channel  leading  to  this  result  is  the  small 
bottom  slope  and  broad,  flat  floodplain  in  which  the  in-bank  channel  meanders. 

Furthermore,  it  was  found  that  the  time  of  arrival  of  the  routed  out¬ 
flow  hydrograph,  in  comparison  with  a  measured  one,  plays  a  significant  role 
in  computing  thereby  T(t);  consequently  computation  of  S(0)  should  not  be 
based  on  the  behavior  of  T(t),  without  taking  the  timing  of  the  routed 
hydrograph  into  account.  In  addition,  the  assumption  that  all  tributary 
inflow  is  concentrated  at  the  lower  end  of  the  reach  can  significantly 


-5- 


affect  the  shape  of  the  desired  T(t),  though  this  factor  is  not  likely  of 
great  importance  in  the  given  reach.  Probably,  the  most  significant  factor 
in  the  deduction  of  negative  tributary  inflow  is  the  assumption  of  an  incor¬ 
rect  storage-outflow  relation  with  too  narrow  a  loop  in  the  Modif ied-Puls 
method . 

In  conclusion,  the  Modif ied-Puls  method  with  a  storage-outflow  based 
on  steady-flow  profiles  is  not  suitable  for  this  reach,  primarily  because 
of  the  great  size  of  the  depth  gradient  relative  to  other  forces  driving  or 
retarding  the  flow,  and  because  the  depth  gradient  is  strongly  a  function 
of  the  unsteadiness  in  a  channel  of  small  slope.  Further,  the  effect  of 
depth  on  storage  is  highly  exaggerated  by  the  broad  flat  floodplains  of 
this  reach.  By  the  same  token,  S(0)  is  highly  event  dependent,  so  that  no 
single  function,  even  approximately  valid  for  a  wide  range  of  peaks,  exists. 

While  doubt  has  been  cast  on  the  accuracy  of  the  topographical  data 
gathered  for  Chuquatonchee  Creek,  it  was  found  that  the  errors  stemming  from 
application  of  Modif ied-Puls  to  events  in  reaches  for  which  it  is  theoretically 
unsuitable,  were  of  sufficient  magnitude  that  they  alone  could  explain  the 
noted  discrepancy. 
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2.  Schema tization  of  the  Chuquatonchee  test  reach  and  flood  events. 

An  examination  of  plots  of  valley  cross  sections  derived  from  large- 
scale  topographical  maps  of  the  Chuquatonchee  tributary  valley  from  Abbott  to 
West  Point,  Mississippi  (drawn  by  the  Mobile  District  of  the  Corps  of  Engineers) 
led  to  schematization  of  the  given  12.4-mile  reach  by  a  prismatic  channel 
characterized  by  the  following  geometrical  parameters. 

In-bank  channel  of  rectangular  section,  width  B  =  70  feet 

depth  D  =  13  feet 

bottom  slope  Sq  =  0.000462  =  2.4  feet/mile 
Manning  n  =  0.07  ft^^ 
length  L  =  12.4  miles 

Single  flood  plain,  characterized  by 

lateral  slope  Spp  =  0.0002 

width  wFp  =  7000  feet 
Manning  n^  =  0.15 

and  vertical  bluffs. 

Historic  flood  events  also  were  schematicized,  two  in  particular:  (1)  the 
small  in-bank  event  of  March  3-9,  1970,  and  (2)  the  large  flood  of  March 
18-24,  1970.  In  both  cases,  inflow-hydrograph  shape  was  assumed  of  Pearson 
type  III,  namely, 

1  1  -  t/tT_ 

I(t)  =  Qq  +  (Ip  -  Qq)  * 

"P 

in  which  Qq  is  the  base  flow,  Ip  is  the  inflow  peak  value,  tip  is  the  time 
to  peak  for  the  inflow  hydrograph,  and  s  is  the  so-called  skew,  held  constant 
in  this  study  at 


s-1 


*  e 


s  -  1 


s  =  1.2 


The  in-bank  event  was  characterized  herein  by  the  peak  value 


I  =  1800  cfs 
1 

with  the  time  to  peak 

t  =48  hours 

P1 

values  similar  to  those  measured  for  event  No.  1. 

The  out-of  banks  flood  was  given  the  values 

I  =  25,000  cfs 
z 

t  =  28  hours 

P2 

again  close  to  the  measured  inflow  characteristics,  for  event  No.  2. 
Base  flow  in  both  cases  was  given  by 
Qg  =  300  cfs. 


The  downstream  boundary  condition  in  both  cases  was  a  normal-depth  rating 
curve . 
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3.  A  preliminary  view  of  Modif ied-Puls  in  the  Chuquatonchee . 

An  approximate  prediction  of  the  efficacy  of  the  Modif ied-Puls  method 
in  the  test  reach  is  made  by  examining  the  dimensionless  comparative  solution 
curves  in  reference  2.  For  a  steady-state  (reference)  base  flow  =  300  cfs, 
the  normal  (reference)  depth  is  Yq=  4.00  feet,  and  the  characteristic 
(reference)  distance  Xq=^q/Sq  =  1.64  miles,  so  that  dimensionless  channel 
length  is  L*  =  L/Xq  =  7.56.  The  reference  time  is  T^  =  xq^qBq/Qq=  2.25  hours. 

For  the  in-bank  event.  No.  1,  with  Ip  =  1800  cfs,  the  ratio  of  peak  to 
base  flow  is  Qp*  =  6.0,  and  the  dimensionless  rise  time  tqp*  =  Tjp  /T^  =  21.3. 
Examination  of  Fig.  lOe  of  ref.  2  suggests  that  Modif ied-Puls ,  with  steps 
fix*  =  4.0  (in  the  given  reach,  with  the  number  of  steps,  N  =  2,  fix*  =  L*/N  = 
3.78)  would  lead  to  outflow  hydrograph  peaks  less  than  2-3%  low  (also,  that 
very  little  attenuation  is  to  be  expected). 

For  the  out-of -banks  event.  No.  2,  Qp*  =  83.3,  and  tqp*  =  12.5.  Exami¬ 
nation  of  Figs.  13e  and  f  of  ref.  2,  for  Qp*  =  5,  suggests  that  Modif ied- 
Puls  outflow  peaks  would  be  about  30%  low.  Comparison  of  Figs.  13d,  for 
Qp*  =  5.0,  and  13g,  for  Qp*  =  40.0,  shows  that  the  percent  error  would 
increase  beyond  the  aforementioned  30%,  for  Qp*  =83.  It  should  be  borne 
in  mind  that  the  curves  of  ref.  2  were  obtained  for  a  base  flow  just  bank 
full,  for  floodplains  with  Manning  n  5  times  that  of  the  main  channel,  and 
for  a  cross  slope  of  0.0019  (0.01  *  Yq/Bq)  .  In  the  given  test  reach,  the 
flood-plain  n  is  just  over  twice  that  of  the  main  channel,  and  the  cross 
slope  is  0.0002.  The  percent-error  figures  must  thus  be  considered  rough 
estimates.  An  indication  of  the  effect  of  the  given  base  flow,  as  opposed 
to  bank  full,  can  be  seen  by  postulating  a  base  flow  of  1800  cfs,  instead 
of  300  cfs.  In  this  case,  Y^=  12.7  feet,  Xq  =  5.2  miles,  T^=  3.78  hours,  so 
that  L*  =  2.38,  Qp*  =  13.9,  tqp*  =  7.4.  With  fix*  =  1,  as  in  Figs.  13d, e 


(as  opposed  to  L*/N  =  1.2',  Modif ied-Puls  peaks  would  be  expected  to  run 
low  by  something  over  20%.  Again,  this  can  only  be  considered  an  estimate 
because  of  the  effects  of  flood-plain  roughness  and  cross  slope. 

An  examination  of  Figs.  1  and  2  shows  that  these  estimates  are  generally 
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Hydrographs  of  event  No.  2  (out-of-banks) 

I  :  Inflow 

0  :  Outflow  (Saint-Venant) 

Og.7  :  Outflow  (Saint-Venant)  with  all  breadths  reduced  30% 

0^  :  Outflow  via  Modif ied-Puls,  N  =  1 

O2  :  Outflow  via  Modif ied-Puls,  N  =  2 

O^g  :  Outflow  via  Modif ied-Puls,  N  =  10 

0^  7  :  Outflow  via  Modif ied-Puls ,  N  =  1,  storage  reduced  to  70% 

O27  •  Outflow  via  Modif ied-Puls ,  N  =  2,  storage  reduced  to  70% 

0?*c;  :  Outflow  via  Modif  ied-Puls ,  N  =  2,  storage  reduced  to  50% 

Outflow  via  Modif ied-Puls  with  base  flow  1800  cfs,  N  =  1 
Outflow  via  Modif ied-Puls  with  base  flow  1800  cfs,  N  =  2 
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4. 


The  role  of  the  storage-outflow  relation  in  the  Modif ied-Puls  method. 


A  typical  storage-outflow  curve  presents  the  former  as  a  monotonically 
increasing  function  of  the  latter,  as  in  Fig.  3. 

s 


0 

Fig.  3.  Typical  storage-outflow  relation. 

With  Eq.  2  in  the  following  form  (tributary  contribution  negligible), 

do  =  L--P  (4) 

dt  dS/dO  v  ' 

and  0  <  dS/dO  <  °°  as  in  Fig.  3,  it  is  clear  that  the  outflow  hydrograph 
must  start  to  rise  immediately ,  without  any  lag,  upon  a  rise  in  inflow. 

This,  in  principle,  is  antithetical  to  the  notion  of  a  wave ,  in  which  an 
event  occurs  at  one  point  in  a  stream,  and  only  some  time  later  is  perceived 
at  a  downstream  point.  Furthermore,  it  is  evident  that  the  peak  outflow 
must  occur  on  the  falling  limb  of  the  inflow  hydrograph,  as  in  Fig.  4. 


Fig.  4.  Relation  between  inflow  and  outflow  hydrographs 
with  "typical"  storage-outflow  relation  (Fig.  3) 


In  actuality,  not  only  does  the  outflow  hydrograph  start  its  rise  some  time 
after  the  inflow  hydrograph,  especially  with  extensive  floodplain  storage, 
but  the  outflow  peak  generally  occurs  somewhat  later  than  at  the  time 
outflow  equals  inflow  (see,  e.g..  Figs.  1  and  2). 

A  delayed  rise  in  outflow  can  occur  (see  Eq.  4)  only  if  dS/dO  is 
infinitely  great  during  the  early  stages  of  flood  inflow.  Similarly,  an 
outflow  peak  at  0  ^  I,  say  0  >  I,  which  is  usually  the  case,  can  also  occur 
only  if  dS/dO  -*•  ®  at  that  time.  Finally,  for  the  outflow  and  inflow  hydro¬ 
graphs  to  cross  without  an  outflow  peak,  i.e.,  0=1,  dO/dt  ^  0,  can  occur 
only  if  dS/dO  =  0  at  that  time.  These  necessary  characteristics  of  the  S(0) 
function  are,  of  course,  precisely  what  are  exhibited  by  the  relations 
computed  by  solution  of  the  Saint-Venant  equations  for  the  schematicized 
Chuquatonchee  reach  and  events  (Figs.  5,  6).1 

In  general,  with  a  given  inflow  hydrograph,  it  is  possible  to  construct 
any  desired  outflow  hydrograph  by  manipulating  the  storage-outflow  relation. 
Some  appreciation  of  the  control  exerted  by  the  S(0)  function  upon  0(t)  can 
be  acquired  by  a  study  of  Fig.  7,  which  shows  in  dimensionless  form  the 
entire  range  of  outflow  hydrographs  that  can  be  obtained  with  a  given  shape 
of  inflow  hydrograph  by  varying  the  slope  of  a  straight-line  S(0)  curve.  All 
discharges  I*,  0*  are  referenced  to  the  inflow  peak  discharge,  and  all  times 
t*  are  referenced  to  the  time  of  peak  inflow.  Dimensionless  storage  S*  is 
consequently  referenced  to  the  product  of  peak  inflow  discharge  and  time  to 

'^'Longitudinal  depth,  area,  and  discharge  profiles  were  obtained  in  the  course 
of  solution  of  these  one-dimensional  hydrodynamic  equations  (by  the  model 
described  in  ref.  2).  This  permits  calculation  of  the  volume  stored  in  the 
reach  at  any  instant.  Coordination  with  the  corresponding  outflow  leads  to 
the  S (0)  function  pertinent  to  the  given  event.  While  S(0)  could  also  have 
been  deduced  by  integrating  dS/dt  =1-0,  the  former  method  illustrates 
quantitatively  how  the  storage  is  distributed  along  the  reach. 


Dimensionless  storage-outflow  relation,  event  No.  1 

a.  Computed  from  a  succession  of  steady  states 

b.  In  effect  during  event  No.  1 

For  storage  in  acre-feet,  multiply  S*  by  55.8;  for  outflow  in  cfs 
multiply  0*  by  300  (ref.  2). 
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peak.  The  dimensionless  slope  of  the  storage-outflow  relation  dS*/dO*  is 
the  theoretical  travel  time  of  the  wave  in  the  reach,  relative  to  the  time 
to  peak  of  the  inflow  hydrograph.  This  value,  t^*,  is  the  parameter  dis¬ 
tinguishing  the  various  outflow  hydrographs  in  Fig.  7.  The  shape  of  the 

inflow  hydrograph  is  again  of  Pearson  type  III  with  skew  s  =  1.2,  and  a 

base  flow  1.2%  of  the  peak.  Relatively  large  increases  in  storage  for  a 
given  increase  in  outflow,  i.e.,  large  values  of  dimensionless  travel  time 
t  *,  lead  to  small,  late  outflow  hydrograph  peaks.  The  smaller  is  tT*,  the 
more  nearly  do  the  inflow  and  outflow  hydrographs  coincide.  In  every  case, 
of  course,  the  area  under  the  complete  hydrograph  is  unity.  That  the  loca¬ 
tions  of  the  outflow  peaks  are  confined  to  the  falling  limb  of  the  inflow 
hydrograph  is  a  consequence  of  the  assumption  of  a  unique  relation  between 
storage  and  outflow,  applied  both  to  the  rising  and  falling  limbs  of  the 

outflow  hydrograph.  Only  a  looped  S(0)  curve,  with  dS/dO  =  0  at  the  time 

that  1=0,  will  allow  the  outflow  hydrograph  to  continue  rising  at  I  =0. 

At  the  time  that  the  outflow  ultimately  peaks,  dS/dO  -+■  <»;  the  subsequent 
descent  of  the  outflow  hydrograph  is  subject  to  that  branch  of  the  loop  that 
descends  monotonically  to  the  starting  value  of  S,  at  the  end  of  the  event. 

As  pointed  out  by  Slocum  and  Dandekar  (ref.  1),  a  single-valued  storage- 
outflow  relation  can  be  made  into  a  loop,  by  applying  the  Modif ied-Puls 
method  as  previously  mentioned,  herein,  to  a  series  of  subreaches  comprising, 
together,  the  given  reach.  The  storage-outflow  relation  deduced  from  the 
resulting  reach-outflow  and  inflow  hydrographs,  say  through  summation  of  the 
increments 

SS  =  (I  -  0)6 1  (5) 

exhibits  a  loop,  the  width  of  the  configuration  depending  upon  the  number  of 


subreaches  N. 


With  the  in-bank  event.  No.  1,  the  relatively  narrow  loop  which  exists 
m  fact  (Fig.  5)  is  fairly  well  modeled  by  the  loop  obtained  with  an  original, 
single-valued  curve  stemming  from  the  assumption  of  normal  depth  for  all 
discharges,  followed  by  application  of  Modified-Puls  to  subreaches,  with  N  =  2. 
The  results  of  the  technique  for  N  =  1  and  2  are  shown  in  Fig.  1. 

For  the  large,  out-of-banks  flood,  event  No.  2,  the  loop  obtained  with 
N  =  2  from  a  similar,  normal-depth,  original  storage-outflow  curve  is 
far  narrower  than  the  true  loop,  shown  in  Fig.  6.  Evidently,  from  the  results 
presented  in  Fig.  2,  N  =  10  results  in  a  loop  of  about  the  right  proportions. 

As  pointed  out  in  ref.  2  it  is  usually  possible  to  find,  empirically, 
a  ''nlue  of  N,  that  will  produce  an  outflow  that  matches  a  given  one  in  peak 
and  time  of  peak.  However,  because  the  storage-outflow  loop  is  produced 
mathematically,  rather  than  physically,  it  is  not  possible  to  predict  the 
optimum  N  from  physical  considerations,  short  of  already  knowing  the  outflow 
hydrograph . 

When  comparing  with  an  actual  outflow  hydrograph  the  results  of  Modified- 
Puls  coupled  with  a  storage-outflow  relation  based  on  area  profiles  derived 
for  a  series  of  steady  flows  and  N  =  1,  the  peaks  of  the  latter  are  seen  to 
be  generally  too  low  and,  of  course,  situated  on  the  descending  limb  of  the 
inflow  hydrograph,  rather  than  somewhat  later.  In  ref.  2,  it  is  shown  that 
increasing  N  both  raises  the  outflow  peak  and  delays  it.  i’ote  how,  in  Fig. 

7.,  the  computed  travel  time  of  the  flow  peak  is  less  than  t^,  a  figure  very 
close  to  the  actual  travel  time  (ref.  2).  This  circumstance,  essentially 
unno ticeable  for  small  travel  times  (t  *  =  0.1,  0.2)  that  result  in  little 
attenuation,  is  greatly  in  evidence  at  travel  times  equal  to  the  time  to 
inflow  peak  and  greater.  Increasing  N,  as  in  ref.  2,  causes  the  computed 
travel  time  of  the  peak,  to  approach  the  theoretical  travel  time  t^,. 


The  computed  magnitude  of  the  peak  outflow  can  also  be  increased  by 
decreasing  the  storage  assumed  at  each  value  of  outflow.  Slocum  and 
Dandekar  (ref.  1)  found  that  uniformly  decreasing  the  storage  by  a  factor 
of  30%  from  the  values  computed  with  HEC-2  for  the  given  reach  of  the 
Chuquatonchee  yielded  about  the  right  values  of  peak  outflow,  computed  by 
Modified  Puls  with  N  =  2.  In  the  schematicized  reach  studied  by  the  writer, 

70%  of  backwater- curve  computed  storage  brought  the  outflow  peak  up  somewhat, 
but  a  factor  of  50%  (with  N  =  2)  was  necessary  to  match  the  correct 
peak  (see  Fig.  2).  Worthy  of  note,  decreasing  the  storage  values  for  a 
given  outflow,  raises  the  peak  outflow,  but  also  speeds  up  the  time  of  arrival 
of  that  peak.  This  plays  a  significant  role  in  the  deduction  of  the  tributary- 
inflow  hydrograph,  as  will  be  discussed  in  section  5. 

It  has  been  suggested  that  the  topographical  data  for  the  Chuquatonchee 
could  be  in  error,  leading  to  HEC-2-computed  volumes  too  great  for  the  given 
outflows.  As  reported  in  ref.  1,  a  70%  factor  was  applied  to  storage  as 
used  in  the  Saint-Venant  routings  as  well.  In  the  present  studies,  in  the 
schematic  channel,  a  30%  reduction  in  width  at  every  depth,  led,  in  the  out- 
of-banks  event,  2,  to  a  3%  increase  in  outflow  peak  computed  with  the  Saint- 
Venant  equations.1  In  the  physically  less  consistent  procedure  applied  to 
the  solution  of  the  Saint-Venant  equations  reported  in  ref.  1,  only  the  top 
width  B  in  the  partial  derivative  of  area  with  respect  to  time,  B9y/9t  = 

9A/9t,  in  the  continuity  equations  was  reduced  by  the  30%.  All  other  terms 
in  the  equations  were  left  undisturbed.  This  procedure,  applied  for  the 
present  study  by  HEC  personnel  to  both  the  schematic  channel  with  event  No.  2 

^When  the  same  width  reduction  was  applied  to  the  in-bank  event,  the  latter 
becomes  an  out-of-bank,  flood  with,  of  course,  totally  different  character- 
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and  the  original  Chuquatonchee  topography  and  corresponding  event,  led  to  a 
10%  increase  in  peak  outflow.  Thus,  the  effect  of  the  topographical  adjust¬ 
ments  suggested  in  ref.  1  are  far  greater  on  the  Modif ied-Puls  results,  than 
on  the  results  of  solution  of  the  Saint-Venant  equations. 

While  the  change  in  Modif ied-Puls  results  due  to  the  70%  storage  factor 
with  real  Chuquatonchee  topography  was  only  about  2/3  of  that  noted  in  the 
schematic  channel,  it  seems  clear  that  this  factor  is  less  a  correction  of 
faulty  topographical  data  than  an  attempt  at  adjusting  a  steady-flow  storage- 
outflow  relationship,  so  it  could  be  used  with  an  unsteady  flow.  Further 
evidence  in  support  of  this  premise  is  found  upon  examination  of  the  Saint- 
Venant-computed  outflow  hydrographs  for  all  five  events  on  the  Chuquatonchee 
reported  in  ref.  1:  most  of  the  peaks  shown,  computed  with  the  70%  factor 
in  force,  run  about  10%  higher  than  the  measured  peaks;  troublesome  negative 
tributary  inflow  is  deduced.  While  other  factors  can  contribute  to  apparent 
negative  tributary  inflow,  it  is  unlikely  that  peaks  correctly  computed  with 
no  tributary  inflow  would  be  higher  than  observed  peaks,  especially  if  this 


inflow  was  in  fact  significant. 


5. 


Deduction  of  tributary  inflow. 


From  the  foregoing  and  Eq.  2  written  in  the  form 


dS  dO 
dO  dt 


T  =  I  -  0 


(6) 


it  is  clear  that  given  inflow  and  outflow  hydrographs  for  any  particular  event 
can  yield  any  T(t)  function  by  choosing  a  commensurate  S(0)  function,  and  vice 
versa.  A  characteristic  feature  of  T(t)  hydrographs  deduced  in  ref.  1  from 
chosen  S(0)  relations  was  the  tendency  of  the  former  to  exhibit  negative 
values  during  certain  portions  of  the  total  time  of  interest.  Several  possible 
causes  of  this  phenomenon  were  investigated. 

The  particular  deduction  technique  employed  in  ref.  1  is  based  on  the 
assumption  that  all  of  the  ungaged  local  flow  is  concentrated  in  a  single 
tributary  entering  the  given  reach  just  above  the  downstream  gage.  Then, 
that  tributary  inflow  hydrograph  is  found  by  subtracting  the  routed  hydrograph 
0R  based  on  zero  local  flows  from  the  observed  outflow  hydrograph  0Q. 

Not  only  the  relative  magnitudes  of  computed  and  measured  outflow  hydro¬ 
graphs  influence  the  inferred  tributary  inflow,  but  also  their  relative 
timing.  Figure  8  shows  the  tributary  hydrographs  deduced  from  three  routed 
hydrographs  0Rl>  0^,  Op^,  differing  only  in  arrival  time,  relative  to  an 
observed  outflow  0Q.  For  simplicity,  the  hydrographs  are  given  triangular 
shapes.  If  Oj^  is  assumed  in  the  correct  position  relative  to  0^,  with  the 
resulting  tributary  hydrograph  T2,  a  routed  hydrograph  which  arrives  too 
early,  such  as  0R^  (or  a  measured  hydrograph  that  appears  to  arrive  later 
than  in  fact)  can  yield  a  tributary  hydrograph  that  is  low  or  negative  in 
its  earlier  portion,  and,  in  compensation,  too  high  in  its  later  portions. 

Of  course  the  net  area  under  each  of  the  deduced  tributary  hydrographs  is 
the  same.  With  a  routed  hydrograph  arriving  later  than  it  should,  relative 


Fig.  8.  Influence  of  outflow  hydrograph  timing  on  deduced  tributary  inflow 
00:  "Observed"  outflow  hydrograph 
0R:  Routed  hydrograph. 

to  the  observed,  excessively  high  values  of  tributary  inflow  are  inferred  at 
small  times,  and  low  or  negative  values,  as  with  Or^,  at  large  times.  Modified 
Puls  routings  tend  to  arrive  too  soon,  so  they  can  be  expected  to  yield 
tributary  inflow  patterns  more  like  Or^  than  Or^.  "Observed"  outflow  hydro¬ 
graphs  tend  to  peak  later  than  in  fact,  because  they  are  generally  based  on 
single-valued  rating  curves.  The  actual  observation  is  of  stage,  which  is 
translated  by  the  rating  curve  to  outflow,  and  in  an  unsteady  flow,  the  stage 
peak  arrives  actually  a  little  later  than  the  discharge  peak.  This  phenomenon. 
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too,  leads  to  tributary-hydrograph  errors  of  the  type,  but  with  overbank 
events  in  the  Chuquatonchee ,  the  magnitude  of  error  from  this  source  is 
relatively  small.  The  width  of  the  loop  in  the  stage-discharge  relation  is 
far  narrower  than  that  of  the  storage-outflow  loop.  Furthermore,  the  errors 
are  somewhat  compensated  in  0  ,  because  a  gaged  inflow  hydrograph  is  subject 
to  the  same  kind  of  error.  Similarly,  the  underestimation  of  0^  (by  a  steady- 
flow  rating  curve)  during  its  rise  in  a  flood,  and  overestimation  during  the 
fall  leads  to  a  relatively  small  T(t)  error,  also  negative  early,  and  positive 
later. 

Topographic  maps  of  the  given  reach  of  the  Chuquatonchee  show  several 
tributaries  entering  the  main  stream  at  various  points  along  the  way.  In 
particular,  a  large  tributary  appears  to  enter  just  downstream  of  the  ap¬ 
parently  smaller,  but  gaged,  Houlka  Creek.  The  contribution  of  the  Houlka 
was  lumped  in  together  with  the  inflow  to  the  upper  end  of  the  Chuquatonchee 
reach  to  comprise  the  inflow  hydrographs  used  in  ref.  1.  Thus,  there  is 
some  concern  about  the  possibility  that  perhaps  most  of  the  tributary  inflow 
occurs  near  the  upstream  end  of  the  reach,  rather  than  at  the  downstream  end 
as  assumed  in  ref.  1  when  deducing  the  tributary-inflow  hydrograph. 

The  influence  of  this  assumption  was  investigated  by  distributing  a 
given  tributary  hydrograph  T(t)  in  various  ways  over  the  length  of  the  reach 
and  noting  the  effect  on  the  outflow  hydrograph.  The  mathematical  model 
used  to  depict  this  behavior  was  againthat  based  on  the  Saint-Venant  equations 
and  described  in  ref.  2.  The  model  was  programmed  for  a  full  range  of 
uniformly  varying  (straight-line)  longitudinal  distributions  of  tributary 
inflow.  Examples  of  these  are  shown  in  Fig.  9.  The  ordinate  scale  therein 


is  the  ratio 


Fig.  9.  Longitudinal  distributions  of  tributary  inflow. 


dT/dx 
r  "  T/L 


(7) 


and  the  parameter  identifying  the  distributions  is  the  upstream  intercept  ry 
or  downstream  intercept  r^.  Of  course,  the  area  under  all  the  lines  is 
unity.  In  addition,  concentrated  inflow  at  either  the  upstream  or  downstream 
ends  (r^  =  °°,  r^  =  ®,  respectively)  was  also  allowed  in  the  model. 

Figure  10a  shows  the  outflow  hydrographs  for  various  distributions  of 
tributary  inflow  augmenting  the  upstream  inflow  of  event  No.  1.  The  instan¬ 
taneous  sum  T  of  tributary  flows  is  given  by  a  Pearson  type  III  distribution 
with  skew  of  1.2,  initial  flow  30  cfs,  10%  that  of  the  channel  base  flow,  a 


Fig.  10a.  Outflow  hydrographs  with  tributary  inflow,  event  No.  1 
0  :  Outflow  with  no  tributary  inflow 

ru:  Upstream  intercept  of  normalized  tributary  distribution 
(see  Fig.  9) 

r^:  Downstream  intercept  of  normalized  tributary  distribution 


t  (hours) 


Fig.  10b.  Deduced  apparent  tributary  inflow,  event  No.  1 
a:  With  r<j  =  8 
b:  With  r,,  =  8. 
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peak  of  9 00  cfs,  3  times  that  of  the  base  flow  and  half  the  upstream- inf low 
peak,  and  a  time  to  peak  of  16  hours,  in  contrast  to  the  48-hour  time  to 
peak  of  the  upstream  inflow  hydrograph.  These  characteristics  of  the  hypo¬ 
thetical  tributary  hydrograph  reflect  those  of  the  tributary  hydrograph 
deduced  in  ref.  1  by  use  of  the  UFP  model  on  the  March  3,  1970  flood.  The 
parameter  distinguishing  one  outflow  hydrograph  from  another  is  the  inter¬ 
cept  of  the  straight-line  distribution,  r^  or  r^.  These  curves  illustrate 
the  effect  of  tributary  distribution  on  the  outflow  hydrograph,  and  can  be 
used  to  show  the  influence  of  this  actual  distribution  on  the  deduced  tribu¬ 
tary  hydrograph,  computed,  as  in  ref.  1,  on  the  assumption  that  all  tributary 
flow  is  concentrated  at  the  downstream  end  of  the  reach.  Subtraction  of  the 
outflow  hydrograph  obtained  without  tributary  inflow,  also  shown  in  Fig.  10a, 
from  the  two  hydrographs  with  the  most  nonuniform  tributary  distributions 
results  in  the  two  curves  of  deduced  tributary  inflow  shown  in  Fig.  10b. 
Evidently,  in  the  translation  of  the  tributary  flow  downstream,  the  peak  is 
moved  slightly  later  in  time  and  higher  in  magnitude. 

These  effects  are  far  more  evident  in  the  overbank  event.  For  this  case, 
the  hypothetical  tributary  hydrograph  was  constructed  similar  to  that  deduced 
in  ref.  1  for  event  No.  2,  the  flood  of  March  18,  1970.  Initial  tributary 
flow  was  20  cfs,  6  2/3%  of  the  base  flow,  the  peak  was  7000  cfs,  28%  of  the 
upstream  inflow  peak,  and  the  time  to  peak  was  18  hours,  in  contrast  to  the 
28  hours  to  peak  of  the  upstream  inflow.  The  tributary  hydrograph  was  again 
given  the  form  of  a  Pearson  type  III  with  skew  of  1.2. 

Figure  11a  shows  hydrographs  of  upstream  inflow  (solid  curve),  inflow 
with  tributary  flow  entering  at  the  upstream  end  (dashed  curve),  outflow 
without  tributary  flow  (solid),  outflow  with  tributary  flow  entering  upstream 
(dashed),  and  outflow  with  tributary  flow  entering  just  above  the  downstream 
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Fig.  11a.  Outflow  hydrographs  with  tributary  inflow,  event  No.  2 
I  :  Inflow  at  upstream  end  of  reach 
I^:  Inflow  at  upstream  end  plus  all  tributary  inflow 
0  :  Outflow  without  tributary  inflow 

0^:  Outflow  with  tributaries  actually  entering  upstream 
0 2-  Outflow  with  tributaries  actually  entering  downstream. 


t  (hours) 


Fig.  lib.  Deduced  apparent  tributary  inflow,  event  No.  2 

Tj :  With  tributaries  actually  entering  upstream 
T2:  With  tributaries  actually  entering  downstream 
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gage  (dot-dash) .  The  resulting  deduced  tributary  hydrographs  (assumed  lumped 
downstream)  are  shown  in  Fig.  lib. 


It  is  evident  that  the  longitudinal  distribution  of  tributary  inflow 
assumed  plays  a  great  role  in  determining  the  shape  of  tributary  hydrograph 
deduced,  but  no  actual,  wholly  positive  distribution  can  lead  to  deduced 
negative  values.  The  source  of  the  latter  must  come  from  other  factors  in¬ 
fluencing  the  routed  hydrograph. 

The  most  important  such  factor  in  the  given  reach  of  the  Chuquatonchee, 
is  the  assumed  storage-outflow  relation.  If  an  inappropriate  relation  is 
used  to  route  the  inflow  hydrograph  downstream,  the  deduced  tributary  hydro¬ 
graph  can  exhibit  all  kinds  of  peculiar  behavior.  To  test  this  premise,  the 
tributary  hydrograph  previously  coupled  to  event  No.  2  was  entered,  in  fact, 
at  the  downstream  end  of  the  channel  and  the  outflow  hydrograph  "observed" 
(i.e.,  computed  by  the  mathematical  model).  This  hydrograph  0Q  subtracted 
from  the  inflow  hydrograph,  yields  the  value  of  time-rate  of  change  of 
storage  minus  tributary  inflow, 

f  -  T  ■  I  -  °0  <8> 

from  Eq.  1.  The  left  side  of  Eq.  8  is  depicted  by  a  solid  curve  on  Fig. 
12a.  Since  the  tributaries  are  in  fact  all  at  the  downstream  end,  they  do 
not  influence  the  storage,  hence  dS/dt  in  Eq.  8  can  be  calculated  by  routing 
the  inflow  hydrograph  without  regard  for  tributary  inflow.  When  this  is 
done  by  the  mathematical  model  to  yield  a  correct  outflow  hydrograph  Or^, 
the  resulting 
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Fig.  12a.  Time  rate  of  change  of  storage  with  and  without  downstream 
tributary  inflow,  event  No.  2 
I  -  Oq:  dS/dt  -  T  (Saint-Venant) 

1 :  dS/dt  (Saint-Venant) 

2:  dS/dt  via  Modif ied-Puls,  N  =  2 

3:  dS/dt  via  modif ied-Puls ,  N  =  2,  storage  reduced  to  70% 
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Fig.  12b. 


Apparent  tributary  inflows  as  functions  of  storage-outflow 
relations,  event  No.  2 

1:  With  correct  storage-outflow  relation  for  event  (Saint-Venant) 
2:  With  normal-depth  storage-outflow  relation,  looped  by  Modified- 
Puls,  N  =  2 

3:  With  70  %  of  normal-depth  storage-outflow  relation,  looped  by 
Modif ied-Puls,  N  =  2. 
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can  be  substituted  in  Eq.  8  to  reproduce  the  correct  tributary-inflow 
hydrograph.  The  left  side  of  Eq.  9  is  plotted  as  a  solid  curve  (1)  in 
Fig.  12a.  The  resulting  deduced  T(t)  is  shown  solid,  (1),  in  Fig.  12b. 

A  second  routed  hydrograph  OR2  was  found  by  the  Modif ied-Puls  method 
using  the  storage-outflow  relation  based  on  steady-flow  profiles  and  N  =  2. 
The  resulting  I  -  Or^  =  dS2/dt  is  plotted  by  dots  and  dashes  (2)  in  Fig. 
12a,  and  the  resulting  T2(t)  by  dots  and  dashes  (2)  in  Fig.  12b.  The 
application  of  the  aforementioned  70%  storage  factor  to  the  Modif ied-Puls 
method  led  to  the  dashed  curves  (3)  of  Figs.  12a  and  b. 


6. 


Conclusions . 


When  the  storage-outflow  relationship  for  use  in  the  Modif ied-Puls 
method  is  computed  in  the  reach  of  Chuquatonchee  Creek  between  Abbott  and 
West  Point,  Mississippi,  by  two  different  methods,  the  results  do  not  agree. 
One  method  computes  the  storage  under  steady-flow  water-surface  profiles  and 
associates  this  with  the  outflow.  The  other  takes  several  past  records  of 
inflow  and  outflow  and  seeks  such  a  storage-outflow  relationship  common  to 
all  these  calibration  events,  that  produces  reasonable  deduced  tributary- 
inflow  hydrographs,  specifically,  one  that  minimizes  their  negative  portions. 

The  investigation  of  this  disagreement  was  performed  through  a  series 
of  easily  controlled  numerical  experiments  in  a  prismatic  schematicized 
channel  similar  in  its  principal  features  to  the  given  Chuquatonchee  reach, 
and  with  schematic  hydrographs  similar  to  the  historic  events.  Real  flows 
in  this  channel  were  mathematically  modeled  by  the  one-dimensional  Saint- 
Venant  equations.  By  this  means,  water-surface  profiles  extant  during  the 
floods  were  computed  and  storage  at  any  instant  readily  determined.  The 
conclusions  drawn  from  the  results  of  the  study  in  the  prismatic  channel 
were  supported  by  computer  runs  using  the  HEC-1,  HEC-2,  and  UFP  programs  with 
real  topographical  and  flood  data  for  the  Chuquatonchee  reach,  as  well  as 
for  the  schematic  reach  and  events.  This  suggests  that  schematization  is  a 
valid  procedure  for  investigating  differences  between  approaches  and  that  the 
nonprismatic  character  of  the  real  channel  does  not  affect  such  results  sig¬ 
nificantly. 

The  principal  reason  for  the  disagreement  between  the  two  methods  of 
determining  a  storage-outflow  relationship  lies  in  the  fact  that  no  such 
(unique)  relationship  exists  for  the  reach.  The  key  postulate  in  the  Modified 
Puls  approach  is  that  there  exists  a  physical  relationship  between  storage 
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and  outflow,  dependent  only  upon  the  topography  and  roughness  of  the  channel. 
If  this  were  true,  the  mathematical  expression  of  this  relation,  S(0),  would 
transcend  a  particular  event,  and  once  found  can  be  used  for  all  events, 
within  the  range  of  calibration.  Instead,  it  was  found  that  the  storage- 
outflow  relation  is  highly  event  dependent. 

The  storage-outflow  relation  computed  through  steady-flow  profiles  is 
in  error,  gross  error  in  the  case  of  overbank  events,  because  the  flow  pro¬ 
files  in  unsteady  flow  are  very  different  from  the  unsteady  ones  at  the  same 
outflow.  The  very  small  bottom  slope  of  the  Chuquatonchee  requires  a  sub¬ 
stantial  depth  gradient  to  move  the  flood  flow.  These  relatively  large  depths 
at  the  upstream  end  of  the  reach  (as  opposed  to  the  small  depths  corresponding 
to  the  outflow  at  the  downstream  end)  engender  an  enormous  amount  of  storage, 
due  to  the  broad  flat  flood  plains  present  in  the  Chuquatonchee  valley. 

Thus,  in  the  early  stages  of  the  flood,  the  storage  in  the  valley  grows  ex¬ 
tensively  while  there  is  practically  no  change  in  outflow.  Then,  eventually, 
when  outflow  (and  downstream  stages)  increase  rapidly  as  the  flood  front 
finally  arrives,  relatively  little  storage  is  added,  and  the  storage-outflow 
relation  is  nearly  flat,'*'  despite  large  changes  in  outflow.  The  trailing 
limb  of  the  outflow  hydrograph  reflects  a  relatively  gradual  decrease  in 
storage  as  the  outflow  falls.  This  results  in  a  very  "obese"  looped  storage- 
outflow  relation.  The  loop  obtained  with  a  single-valued  S(0)  and  N  =  2  is 
not  nearly  fat  enough  to  model  overbank  events.  Further,  while  Modif ied-Puls 
outflow  peaks  can  be  increased  by  arbitrarily  reducing  the  steady-flow- 
computed  storage  for  every  outflow,  the  loop  is  still  too  narrow  and  leads 

^Storage  plotted  as  ordinate,  outflow  as  abscissa. 
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to  peculiar  behavior  of  deduced  tributary  inflow  hydrographs,  including 
double-humps,  regions  of  negative  inflow,  and  so  on. 

Storage-outflow  relationships  computed  from  observed  inflow  and  outflow 
hydrographs  also  do  not  serve  the  purpose,  in  principle,  because  each  dif¬ 
ferent  event  is  characterized  by  a  different  relation,  and  in  application, 
because  of  the  single-valuedness  or  narrow  loop  assumed  for  the  sought-after 
relation.  Any  desired  tributary-hydrograph  behavior  can  be  achieved  with 
given  inflow  and  outflow  hydrographs,  by  adjusting  the  S(0)  relation 
up  to  the  point  of  maximum  outflow.  Then  the  die,  so-to-speak,  has  been 
cast,  and  the  trailing  limb  of  the  routed  hydrograph  is  "at  the  mercy"  of 
the  previously  set  S(0),  with  attendant  loss  of  control  over  that  portion  of 
the  tributary  hydrograph. 

Concern  over  the  accuracy  of  the  topographical  data  for  the  Chuquatonchee 
became  secondary,  when  it  was  noted  that  these  theoretical  considerations 
were  sufficient  in  themselves  to  produce  the  difficulties  encountered  in 
routing  floods  in  this  reach.  When  the  geometrical  data  was  adjusted  to 
decrease  storage  by  30%,  in  an  earlier  study,  then  used  in  the  UFP  program, 
the  resultant  outflow  peaks  computed  exceeded  the  measured  values,  by  about 
10%.  The  tests  performed  for  the  present  study  suggest  that  use  of  the 
original  topographic  data  might  well  reduce  these  to  observed  levels  and 
reduce  the  incidence  of  deduced  negative  tributary  inflow. 

Similarly,  errors  in  estimation  of  roughness  would  be  relatively  insig¬ 
nificant.  Other  factors,  such  as  longitudinal  distribution  of  tributary 
inflow,  and  errors  in  inflow  and  outflow  measurements  due  to  unsteadiness 
were  investigated  and  found  to  be  relatively  unimportant,  compared  to  the 
clear  event  dependency  of  the  storage-outflow  relation. 
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